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ABSTRACT 

We investigate the effect of radiative feedback on the star formation process using radia- 
tion hydrodynamical simulations. We repeat the previous hydrodynamical star cluster forma- 
tion simulations of Bate et al., and Bate & Bonnell, but we use a realistic gas equation of 
state and radiative transfer in the flux-limited diffusion approximation rather than the original 
barotropic equation of state. 

Whereas star formation in the barotropic simulations continued unabated until the simula- 
tions stopped, we find that radiative feedback, even from low-mass stars, essentially terminates 
the production of new objects within low-mass dense molecular cloud cores after roughly 
one local dynamical time. Radiative feedback also dramatically decreases the propensity of 
massive circumstellar discs to fragment and inhibits fragmentation of other dense gas (e.g. 
filaments) close to existing protostellar objects. These two effects decrease the numbers of 
protostars formed by a factor of « 4 compared with the original hydrodynamical simulations 
using the barotropic equation of state. In particular, whereas the original simulations produced 
more brown dwarfs than stars, the radiative feedback results in a ratio of stars to brown dwarfs 
of approximately 5:1, in much better agreement with observations. Most importantly, we find 
that although the characteristic stellar mass in the original calculations scaled linearly with the 
initial mean Jeans mass in the clouds, when radiative feedback is included the characteristic 
stellar mass is indistinguishable for the two calculations, regardless of the initial Jeans mass 
of the clouds. We thus propose that the reason the observed initial mass function appears to 
be universal in the local Universe is due to self-regulation of the star formation process by 
radiative feedback. We present an analytic argument showing how a characteristic mass may 
be derived that is relatively independent of initial conditions such as the cloud's density. 

Key words: stars: formation, stars: low-mass, brown dwarfs, stars: luminosity function, mass 
function, methods: numerical, radiative transfer, hydrodynamics. 



1 INTRODUCTION 

Understanding the origin of the stellar initial mass function (IMF) 
is one of the fundamental goals of a complete theory of star for- 
mation. Observationally, much has been accomplished in the five 
decades since Salpeter ( 1955) published his seminal paper on the 
form of the IMF. Salpeter's determination of the high-mass end 
of the IMF has become widely accepted, but the behaviour of the 
IMF below 1 Mq has been more accurately characterised, with the 
turnover into the brown dwarf regime being the subject of many 
recent investigations. The general form of the IMF in the solar 
neighbourhood is now known down to ~ 0.03 Mq ( Krou pa[2001| 
Chabrier 2003 ). However, despite this progress in determining the 
form of the IMF, there is still no standard model for the origin of 
the IMF or on how it should depend on environment. In fact, much 
of the difficulty in determining which processes are responsible for 
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the origin of the IMF can be attributed to the fact that the IMF 
doesn't seem to vary strongly with environment (for recent reviews 
see Kroupa 2002 ; Chabrier 2003 ). Galactic studies and those ex- 
tending to the Magellanic Clouds have repeatedly failed to find any 
systematically robust and statistically significant variation from the 
general form of the IMF. There is evidence from observations of 
the Arches Cluster and the stars orbiting the supermassive black 
hole, Sgr A*, that the IMF near the Galactic Centre may be biased 
in favour of massive stars. However, the apparently top-heavy IMF 
in the Arches Cluster jFiger et al-1|1999[ |Stolte et al?]|2002| > may 
be due to dynamical evolution rather than being primordial in ori- 



gin (Portegi es Zwart et al.|2002||Kim et al.|2006||Portegies Zwart| 
|et al.|2007 1, and the apparent deficit of low-mass stars surround- 
ing Sgr A* is only indirectly inferred from the lack of X-ray emis- 
sion ( |Nayakshin" & Sunyaev 2005). This has led to the use of the 
term 'universal IMF' when describing present-day star formation. 
Only for the first, zero metallicity, stars in the Universe is there gen- 
eral agreement that the IMF should differ from that found locally, 
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producing substantially more massive stars (e.g., |Bromm, Coppi &| 
|Larson|199"9||Abel, Bryan & Norman|2 000). However, this has yet 
to be confirmed observationally. 

Many theories have been proposed for the origin of the IMF. 
These fall into four main classes. The IMF may originate from 
fragmentation, whether it be turbulent fragmentation (Larso n[T992| 
Henriksen|[T986l [T99T] |Elmegreen| [T997] [T999] |2000[ IPadoanll 
Nordlund & Jones|1997{|Padoan & Nordlurid|2002}, g ravitational 



fragmentatio n (|Larson 1973| [Elmegreen & Mathieu||1983| |Zin- 
necker]|T984j |Yoshii & Saio||1985| >, or domain packing ( |Richtler| 



1994[>, with the fragmentation subject to an opacity limit which sets 



a minimum stellar mass ( Hoyl e1[T953| |Gaustad||1963[ |Yon eyama 
Wf2\ |Suchkov & Shchekmov||1976||Low & Lynden-Bell||1976| 
Rees|l976||Silk|1977a|b[|Masunaga & Inutsuka|1999| ). It may de- 



pend on feedback processes such as winds and outflows ( [Shu et ak] 
[T9881 |Silk|1995[ | Adams & Fatuzzo|1996J or heatin g of accretion 
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[T9781 |Zinnecker 1982, Bon nell et al.|1997||200T 
[200Tp ; [Kiessen, Burkert & Bate|1998"l|Myers|2000). Or it may be 
due to coalesence or collisional build-up dSilk & Takahash ijl979| 
Pumphrey & Scalo|1983] |Bonnell, Bate & Zinnecker|1998[|Bon 



nell & Bate 2002). In reality, all of these processes are likely to 
play some role. The main questions to answer are which process, if 
any, dominates the origin the the IMF, and how does the IMF vary 
with environment? 

The most fundamental parameter of the IMF is its character- 
istic mass, ~ O.3M0. A lower limit to the mass of a 'star' is pro- 
vided by the opacity limit for fragmentation of a few Jupiter masses 
for typical molecular gas (Hoyle 1953, |Low & L ynden-Bel l|l976| 
|Rees| 1 976[ [Whitworth & Stamatellos 2006). All objects must have 
masses greater than this minimum mass. However, the peak of the 
IMF occurs at masses roughly two orders of magnitude greater than 
this minimum mass. Several theories for the IMF link the charac- 
teristic mass to the typical Jeans mass in the progenitor molecular 
cloud. This may be the thermal Jeans mass (e.g. |Larson|19 92 ), a 
magnetic critical mass, or a turbulent Jeans mass ( Silk 1995J). A 
Jeans mass origin for the characteristic stellar mass has been sup- 
ported by some hydrodynamical calculations of the fragmentation 
of clumpy and turbulent molecular clouds. In these calculations, it 
was found that the mean mass of the protostars was similar to the 
initial mean Jeans mass of the cloud dKlessen et al. 1998 , KlesseiT&l 
|Burkert|2000l[2W][Kie"ssen|2001[[B"ate, Bonnell & Bromm|2003) 
and that variations in the initial Jeans mass led to corresponding 
variations in the characteristic mass of the IMF (Bate & Bonnell 
[20051 [Bate] 2005 •• Another quite different model of the IMF pro- 
poses that the IMF orginates from the mass distribution of dense 
cores in turbulent molecular clouds ( Padoan & Nordlund 2002 1. In 
this model, the characteristic mass of the IMF depends on both the 
initial Jeans mass and the Mach number of the turbulence. How- 
ever, as pointed out by Adams & Fatuzzo 1 1996|>, t here is no unique 
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Jeans mass in a molecular cloud, 
that the appropriate Jeans mass may be that at which gas and grains 
couple thermally and dust cooling takes over from atomic line cool- 
ing in molecular clouds. Whitworth, Boffin & Francis ( 1998} also 
link the characteristic mass of the IMF with this gas-grain coupling 
point. Hydrodynamical simulations using an equation of state in- 
spired by Larson's models and varying it to give a different charac- 
teristic Jeans mass do produce corresponding changes in the char- 
acteristic mass of the IMF ( Jappsen et al. 2005 ; Bonnell, Clarke & 
|Bate|20 06). Recently, Elmegree n, Klessen & Wilsonj ( |2008| have 
argued that the characteristic Jeans mass at which the gas-dust cou- 



pling occurs is relatively independent of environmental quantities 
such as density, temperature, metallicity, and the radiation field, 
perhaps explaining the apparent universality of the characteristic 
mass of the IMF. 

In this paper, we investigate the effect of radiative feedback 
on the star formation process. We repeat the hydrodynamical cal- 
culations of star cluster formation Bate et al. (2003) and |Bate~&| 
Bonnell (2005) that resolved fragmentation down to the opacity 
limit, circumstellar discs, and binary stars (hereafter BBB2003 and 
BB2005, respectively). The initial conditions for these two calcu- 
lations were identical except that the latter cloud was denser than 
the former, lowering the initial mean Jeans mass by a factor of 3. 
BB2005 showed that the characteristic masses of the two IMFs pro- 
duced depended linearly on the initial Jeans mass, also differing by 
a factor of 3. Here, we repeat the original calculations, but instead 
of using a barotropic equation of state to model the thermodynamic 
behaviour of the gas, we use a realistic gas equation of state and in- 
clude radative transfer in the flux-limited diffusion approximation. 
We find that the inclusion of radiative feedback from the forming 
protostars substantially weakens the dependence of the character- 
istic mass of the IMF on the initial Jeans mass in the progenitor 
molecular cloud. Thus, we propose that star formation regulates it- 
self via radiative feedback to provide the IMF with a characteristic- 
mass that is usually only weakly dependent on environment. 

The paper is structured as follows. In Section[2] we briefly de- 
scribe the numerical method and the initial conditions for the sim- 
ulations. We present the results from our radiation hydrodynamical 
simulations in Section[5] In Section]?] we compare our results with 
previous simulations and with observations, and we discuss the im- 
plications of our results for the origin of the characteristic mass of 
the IMF. Our conclusions are given in Section[5] 



2 COMPUTATIONAL METHOD 

The calculations presented here were performed using a three- 
dimensional SPH code based on a version originally developed 
by Benz ( |Benz|1990"[|Benz et al.|1990|> but substantially modified 
as described in Bate et al. ( 1995), Whitehouse, Bate & Monaghan 
P005^ and |Whitehouse & BateH2006) and parallelised by M. Bate 
using OpenMP. 

Gravitational forces between particles and a particle's near- 
est neighbours are calculated using a binary tree. The smoothing 
lengths of particles are variable in time and space, subject to the 
constraint that the number of neighbours for each particle must re- 
main approximately constant at iVneigh = 50. The SPH equations 
are integrated using a second-order Runge-Kutta-Fehlberg integra- 
tor with individual time steps for each particle (jBate et al.|1995 i. 
We use the standard form of artificial viscosity ( Monaghan & Gin- 
|gold| 1 983 ; Monaghan 1992 ) with strength parameters a v = 1 and 
/?v = 2. 



2.1 Equation of state and radiative transfer 

The original hydrodynamical calculations of BBB2003 and 
BB2005 modelled the thermal behaviour of the gas without per- 
forming radiative transfer using a barotropic equation of state. 

In this paper, we use an ideal gas equation of state p — 
pTlZ/jj., where p is the density, T is the gas temperature, 1Z is the 
gas constant, and p is the mean molecular mass. The equation of 
state takes into account the translational, rotational, and vibrational 
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degrees of freedom of molecular hydrogen (assuming an equilib- 
rium mix of para- and ortho-hydrogen; see Black & Bode nheimer| 
|1975| >. It also includes the dissociation of molecular hydrogen, and 
the ionisations of hydrogen and helium. The hydrogen and helium 
mass fractions are X — 0.70 and Y — 0.28, respectively. The 
contribution of metals to the equation of state is neglected. Further 
details on the implementation of the equation of state can be found 
in |Whitehouse & Bate| < [2006l >. 

Two temperature (gas and radiation) radiative transfer in the 
flux-limited diffusion approximation is implemented as described 
by |Whitehouse et al.|p005) and |Whitehouse & Bate| < p006) . En- 
ergy is generated when work is done on the gas or radiation fields, 
radiation is transported via flux-limited diffusion, and energy is 
transferred between the gas and radiation fields depending on their 
relative temperatures, the gas density, and the gas opacity. We use 
interpolation from the opacity tables of Poll ack et al.| (T985) to pro- 
vide the interstellar grain opacities for solar metallicity molecular 
gas, whilst at higher temperatures when the grains have been de- 
stroyed we use the tables of Alexander ( 1975 ) (the IVa King model) 
to provide the gas opacities (for fur ther details, see|Whitehouse &| 
|Bate|2006) . The only change from |Whitehouse & Bate| ( |2006} is 
to do with the boundary condition. The clouds modelled here have 
free boundaries. To provide a boundary condition for the radiative 
transfer, all particles with densities less than 10~ 21 g cm~ 3 have 
their gas and radiation temperatures set to the initial values of 10 
K. This gas is 2 — 3 orders of magnitude less dense that the initial 
clouds (see below) and, thus, these boundary particles surround the 
region of interest in which the star cluster forms. 

2.2 Sink particles 

In the original hydrodynamical calculations of BBB2003 
and BB2005, gas collapsed isofhermally until a density of 
10 -13 g cm -3 beyond which point the gas temperature increased 
with density as p 2 ^ 5 . This resulted in the formation of pressure- 
supported fragments with initial masses of a few Jupiter-masses 
in collapsing regions. As each fragment accreted its central den- 
sity and temperature increased, resulting in smaller and smaller 
timesteps. When the central density of a fragment exceeded 
10 -11 g cm , it was replaced by a sink particle. 

In this paper, the evolution a collapsing region of gas is sim- 
ilar except that the temperature is calculated self-consistently us- 
ing radiative transfer and the realistic equation of state and we fol- 
low fragments to higher densities before inserting a sink particle. 
Collapsing gas that exceeds ~ 10 -13 g cm -3 becomes optically 
thick and traps radiation, thus heating up and forming a pressure- 
supported fragment known as the 'first hydrostatic core' ( |Larson| 
|1969| >. This core continues to accrete mass and its central density 
and temperature rise until molecular hydrogen begins to dissociate 
at a temperature of w 2000 K and density of ~ 10~ 7 g cm -3 . 
The dissociation of molecular hydrogen requires energy, initiating 
a second dynamical collapse within the first core. When the density 
exceeds ~ 10 -3 g cm~ 3 the gas is atomic and a second, 'stellar', 
core is formed. We follow the gas through the entire first core phase 
and through most of the second collapse phase, inserting a sink par- 
ticle when the density exceeds 10 -5 g cm -3 (in terms of the real 
star formation process, this is just a few days before the stellar core 
is formed). 

In the original barotropic calculations, a sink particle was 
formed by replacing the SPH gas particles contained within r acc = 
5 AU of the densest gas particle by a point mass with the same 
mass and momentum. Any gas that later fell within this radius was 



accreted by the point mass if it was bound and its specific angular 
momentum was less than that required to form a circular orbit at 
radius r acc from the sink particle. Thus, gaseous discs around sink 
particles could only be resolved if they had radii J> 10 AU. Sink 
particles interacted with the gas only via gravity and accretion. The 
angular momentum accreted by a sink particle was recorded but 
played no further role in the calculation. The gravitational accel- 
eration between two sink particles was Newtonian for r 4 AU, 
but was softened within this radius using spline softening ( |Benz| 
1990). The maximum acceleration occurred at a distance of « 1 
AU; therefore, this was the minimum separation that a binary could 
have even if, in reality, the binary's orbit would have been hard- 
ened. 

In this paper, our default is to use accretion radii of only 
face = 0.5 AU and we do not use any gravitational softening be- 
tween two sink particles. This allows us to resolve smaller discs 
and closer binaries than in the earlier hydrodynamical simulations. 
We also performed a calculation based on the BBB2003 initial con- 
ditions with accretion radii of 7- acc = 5 AU to see what difference 
this made to the results. Sink particles are permitted to merge if 
they pass within 0.02 AU of each other (i.e., w 4 Rq). This radius 
was chosen because recently formed protostars are thought to have 
relatively large radii (e.g., Larson 1969). No mergers occurred in 
any of the calculations reported here. 

Finally, we emphasise that sink particles only interact with the 
rest of the simulation via gravity and by accreting gas. In partic- 
ular, for the radiative hydrodynamical calculations presented here 
there is no radiative feedback from the sink particles into the rest 
of the simulation. Neglecting the intrinsic luminosity of an accret- 
ing low-mass protostar is reasonable since the accretion luminosity 
overwhelms the intrinsic luminosity of the object even for a very 
low accretion rate. However, since gas is not modelled within r acc 
of each protostar due to the sink particle approximation, a substan- 
tial fraction of the accretion luminosity is also neglected (from the 
gravitational energy liberated during the inspiral of gas from r aC c 
to the stellar surface). Thus, we emphasise that these calculations 
give a lower limit on the effects of radiative feedback on the star 
formation process. We perform one of the calculations with accre- 
tion radii of both 5 AU and 0.5 AU specifically to investigate the 
effect of this approximation on the results. 

2.3 Initial conditions 

The initial conditions for the two calculations are identical to those 
of BBB2003 and BB2005 and are summarised in Table [T] Each 
spherical cloud contains 50 Mq of molecular gas. The radii of the 
two clouds are 0.188 pc and 0.090 pc, respectively, so that the latter 
cloud has 9 times the density of the former (densities of 1.2 x 10~ 19 
and 1.1 x 10~ 18 g cm -3 , respectively). At the initial temperature 
of 10 K, the two clouds have mean thermal Jeans masses of 1 and 
1/3 Mq, respectively. Although each cloud is uniform in density 
initially, an initial supersonic 'turbulent' velocity field is imposed 
on each of them in the same manner as Ostriker, Stone & Gammie 
(2001). The divergence-free random Gaussian velocity field has a 
power spectrum P(k) oc fc -4 , where k is the wavenumber. In three 
dimensions, this results in a velocity dispersion that varies with dis- 
tance, A, as <x(A) oc A 1 / 2 in agreement with the observed Larson 
scaling relations for molecular clouds (Larson 1981 1. The velocity 
field was generated on a 128 3 uniform grid and the velocities of the 
particles were interpolated from the grid. The same velocity field is 
used for each of the two clouds, but the normalisation differs and is 
set so that the the kinetic energy of the turbulence equals the mag- 
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Table 1. The initial conditions and the statistical properties of the stars and brown dwarfs formed in the original BBB2003 and BB2005 calculations and 
the new versions of those calculations using radiation hydrodynamics that we present here. Calculations BBB2003 RT0.5 and BB2005 RT0.5 use radiation 
hydrodynamics and sink particles with accretion radii of r acc = 0.5 AU. Calculation BBB2003 RT5 is identical to BBB2003 RT0.5 except that the sink 
particles have r acc = 5 AU. The initial conditions for the two different types of calculation are identical except that the BB2005-type initial cloud has a 
smaller radius giving a density 9 times higher and a mean thermal Jeans mass a factor of 3 lower. All calculations were run for 1 .40 initial cloud free-fall 
times. Brown dwarfs are defined as having final masses less than 0.075 Mq . The numbers of stars (brown dwarfs) are lower (upper) limits because some of the 
brown dwarfs were still accreting when the calculations were stopped. Using radiation hydrodynamics dramatically reduces the numbers of objects formed, 
particularly brown dwarfs. Furthermore, whereas using a barotropic equation of state led to the median mass scaling linearly with the mean thermal Jeans mass 
of the cloud, radiation hydrodynamics removes the dependence of the median stellar mass on the initial Jeans mass. 



nitude of the gravitational potential energy of each cloud. Thus, the 
initial root-mean- square (rms) Mach number of the turbulence was 
M = 6.4 in BBB2003 and M = 9.2 in BB2005. The initial free- 
fall times of the two clouds are ts = 6.0 x 10 12 s or 1.90 x 10 5 
years and tg — 2.0 x 10 12 s or 6.34 x 10 4 years, respectively. 

In fact, the early evolution of each of the clouds was not re- 
peated using radiation hydrodynamics since the gas remains essen- 
tially isothermal during the early evolution with radiative transfer, 
but the radiation hydrodynamical calculations are approximately 
an order of magnitude more computationally expensive than the 
barotropic calculations. The radiation hydrodynamical calculations 
were instead begun from dump files from the original calculations 
just before the density exceeded 1 x 10 -16 g cm -3 and 9 x 10~~ 16 
g cm~ 3 , respectively. 



2.4 Resolution 

The local Jeans mass must be resolved throughout the calculations 
to model fragmentation correctly (Bate & Burkert 1997| |Truelove 
et al.|1997|rWhitworth|1998]|Boss et al |2000| |Hubber, G oodwin 



& Whitworth|2006| >. This requires J> 1.5iV nc jgh SPH particles per 
Jeans mass; jV ne i g h is insufficient (BBB2003). The original calcu- 
lations used 3.5 x 10 6 particles to model the 5O-M.0 clouds and 
resolve the Jeans mass down to its minimum value of ~ 0.0011 
Mq (1.1 Mj) at the maximum density during the isothermal phase 
of the collapse, p cr it = 10~ 13 g cm -3 . Using radiation hydro- 
dynamics results in temperatures at a given density no less than 
those given by the original barotropic equation of state and, thus, 
the Jeans mass is also resolved in the calculations presented here. 

Each of the two calculations with sink particle accretion radii 
of r acc = 0.5 AU required the equivalent of approximately 40,000 
CPU hours on a 16-processor 1 .65GHz IBM p570 compute node of 
the United Kingdom Astrophysical Fluids Facility (UKAFF). One 
of the two calculations was completed on the University of Exeter 
Supercomputer, an SGI Altix ICE 8200. The r acc = 5 AU calcula- 
tion ran approximately 20 times faster. 



3 RESULTS 

Results from three new calculations are presented in this pa- 
per. All were performed using the radiation hydrodynamical SPH 



code described above. Calculations BBB2003 RT0.5 and BB2005 
RT0.5 are identical the hydrodynamical calculations presented in 
BBB2003 and BB2005, respectively, except that they were per- 
formed using radiation hydrodynamics and sink particles have ac- 
cretion radii of 0.5 AU, ten times smaller than in the original calcu- 
lations. Calculation BBB2003 RT5 is identical to BBB2003 RT0.5 
except that the sink particles have accretion radii of 5 AU, the same 
size as in the original BBB2003 calculation. This last calculation 
enables us to determine how sensitive the results are to the size 
of the accretion radii. This is even more important when using ra- 
diation hydrodynamics than for a barotropic equation of state be- 
cause sink particles themselves do not emit radiation - only the 
gas emits radiation (see Section [2~2| . Thus, using smaller accretion 
radii allows more accretion luminosity to be released by the pro- 
tostars into the calculations, leading to hotter gas and potentially 
affecting the pattern of fragmentation. Each of the radiation hy- 
drodynamical simulations is followed to 1.40 iff, the same as the 
original BBB2003 and BB2005 calculations. 



3.1 BBB2003 initial conditions 

We begin by presenting the results from the two BBB2003-type 
calculations using radiation hydrodynamics. As mentioned in Sec- 
tion |2.3| the radiation hydrodynamical calculations were not re-run 
from the initial conditions, but were started from the last dump 
file from the original BBB2003 before the density exceeded 10 -16 
g cm -3 . Before this point the initial 'turbulent' velocity field had 
generated density structure in the gas, some of which was collected 
into dense cores which had begun to collapse. Those readers inter- 
ested in this early phase should refer to BBB2003 for figures and 
further details. 

The BBB2003 radiation hydrodynamical calculations were 
started from t = 0.976 ts (in the original BBB2003 calculation 
the first sink particle was inserted at t — 1.037 ts, some 1.2 x 10 4 
years later). Using radiation hydrodynamics, the first sink particle 
is inserted at t — 1.042 tg. The slightly later time is primarily be- 
cause in the original calculation sink particles were inserted when 
the density exceeded 10~ 41 g cm~ 3 (when the fragment was in the 
'first hydrostatic core' stage of evolution) whereas with the radi- 
ation hydrodynamics we do not insert sink particles until halfway 
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Figure 1. Global evolution of the BBB2003 RT0.5 (left) and BB2005 RT0.5 (right) calculations. In each case, the left (red-yellow-white) panels show the 
logarithm of column density, N, in g cm -2 and the scales cover — 1.7 < log N < 0.5 (left) and —1.1 < log N < 2.1 (right). The right (blue-red-yellow- 
white) panels show the logarithm of mass weighted temperature, T, in K with the scales covering 9 — 100 K. The lefthand panels are 0.4 pc (82,400 AU) 
across while the righthand panels are 0.192 pc (39,600 AU) across. Time is given in units of the initial free-fall time of 1.90 X 10 5 yr (left) and 6.34 X 10 4 
vr (right). 
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Figure 2. The star formation in the main dense core of the BBB2003 RT0.5 calculation. The first object forms at t = 1.042 tg. Large gaseous filaments 
collapse to form single objects and multiple systems. These objects fall together to form a small group. Radiative feedback from the accreting protostars heats 
the gas in the dense core. Each panel is 0.025 pc (5150 au) across. These may be compared to the equivalent figures in the original BBB2003 paper. Time is 
given in units of the initial free-fall time of 1.90 X 10 yr. The red-yellow- white panels show the logarithm of column density, N, through the cloud, with the 
scale covering —0.5 < log TV < 2.5 with N measured in g cm -2 . The blue-red-yellow-white panels show the logarithm of mass weighted temperature, T, 
through the cloud with the scale covering 9 — 300 K. 



through the second collapse phase at a density of 10 5 g cm 3 (see 
|2.2|for further details). 



3.1.1 Sink particles with r acc = 0.5 All 

In the lefthand panels of Figure [T] we present snapshots of the 
global evolution of the BBB2003 RT0.5 calculation. The left 
most panels (using the red-yellow-white colour scale) display 
the column density, while the centre-left panels (using the blue- 
red-yellow-white colour scale) display the mass weight temper- 
ature in the cloud. Figure [2] shows the evolution of the main 
dense core from the BBB2003 RT0.5 calculation in much greater 



detail than Figure ^ These snapshots are shown at the same 
times as the equivalent figures in the original BBB2003 pa- 
per. An animation comparing the original calculation with the 
radiation hydrodynamical calculation can be downloaded from 
http://www.astro.ex.ac.uk/people/mbate~ 

Comparison of the snapshots and/or the animations shows that 
the barotropic and radiation hydrodynamical calculations diverge 
quickly on small scales. In the original calculation, the first pro- 
tostar to form is surrounded by a massive circumstellar disc that 
quickly fragments into three more objects - two brown dwarfs and 
a low-mass star. With radiation hydrodynamics, this massive disc 
does not fragment. The accretion luminosity released as gas falls 
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Figure 3. The star formation in the second (left) and third (right) dense cores at the end of the BBB2003 RT0.5 calculation. The left two panels are 600 AU 
across, while the right two panels are 1000 AU across. These may be compared to the equivalent figures in the original BBB2003 paper. Time is given in units 
of the initial free-fall time of 1.90 X 10 5 yr. The red-yellow-white panels show the logarithm of column density, N, through the cloud, with the scale covering 
0.0 < log N < 2.5 with N measured in g cm -2 . The blue-red-yellow-white panels show the logarithm of mass weighted temperature, T, through the cloud 
with the scale covering 9 — 300 K. At the end of the calculation, the object in the top-right corner of the panels depicting the third core had not yet undergone 
the second collapse and been replaced by a sink particle. 
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Figure 4. The cumulative distributions of the distance between a new pro- 
tostar and its closest other protostar (excluding the first protostar to form in 
each calculation) for each of the 5 calculations discussed in this paper. The 
distributions from the previously published BBB2003 and BB2005 calcula- 
tions using a barotropic equation of state are given by the thin solid line and 
thin dashed line, respectively. The radiation hydrodynamical calculations 
presented here are BBB2003 RT0.5 (thick solid line), BB2005 RT0.5 (thick 
dashed line), and the large accretion radius BBB2003 RT5 calculation (thin 
dot-dashed line). It is clear that including radiative feedback increases the 
typical distance between a new protostar and its closest companion. For the 
BB2005-type initial conditions, 3/4 of the objects formed within 300 AU of 
another protostar using a barotropic equation of state, whereas with radia- 
tive feedback this fraction is reduced to less than 1/3. For the BBB2003-type 
initial conditions, more than 70% of the objects formed within 500 AU of 
another protostar using a barotropic equation of state, whereas with radia- 
tive feedback and small accretion radii this percentage is reduced to just 
over 40%. 



onto the disc and then spirals in towards the central protostar is 
sufficient to heat the disc and prevent its fragmentation. This is 
one of the two main differences between the original calculation 
using the barotropic equation of state and the radiation hydrody- 
namical calculation - discs and dense filaments close to existing 
protostars are inhibited from fragmenting by the radiative feedback 
due to the accretion luminosity released by the low-mass stars and 
brown dwarfs. This is not surprising. Whitehouse & Bate (20061 



previously showed that replacing the barotropic equation of state 
with radiative transfer can lead to temperatures up to an order of 
magnitude higher near young low-mass protostars and, thus, po- 
tentially strongly inhibits fragmentation, while Krumholz (2006 1 



made a similar argument analytically. Rafikov (2005 1 



Matzner & 



Levml p005l >, |Kratter & Matzner| j2006| >, and [Whitworth &"Sta 



matellos (20061 have all pointed out that accurate treatments of 
radiative transfer are likely to significantly decrease disc fragmen- 
tation. Finally, we note that the temperature field surrounding the 
newly-formed protostars varies significantly on timescales of hun- 
dreds to thousands of years. In the animations, this can be seen as 
'flicking' of the temperature field. These temperature variations are 
due to variations in the accretion rates of the protostars and their 
discs, particularly when protostars undergo dynamical interactions 
that perturb their discs. 

The radiative feedback primarily affects fragmentation on 
length scales ^ 500 AU. As can be seen by comparing the snap- 
shots in Figure[2]with the equivalent figures in BBB2003, the larger 
scales are unaffected, at least initially. In Figure [4] we plot the cu- 
mulative distribution of the distances between each object and its 
closest other protostar at the time of formation of the object (i.e., 
the time a sink particle is inserted). The thin solid line gives the 
distribution for the original barotropic BBB2003 calculation, while 
the thick solid line gives the distribution for the BBB2003 RT0.5 
calculation. Comparison of the two distributions shows that the 
fraction of objects forming within 100 AU of an existing proto- 
star is twice as large without radiative feedback. Similarly, more 
than 80% of the objects form within 1000 AU of an existing pro- 
tostar in the barotropic calculation while this percentage is reduced 
to 50% with radiative feedback. The small-scale effect of radiative 
feedback is in contrast to the effects of magnetic fields (see |Price &| 
|Bate|20 08 , who performed simulations similar to BBB2003 with a 
range of magnetic field strengths) which affect the cloud's structure 
on both large and small scales. Eventually, however, radiative feed- 
back begins to alter structures on large scales indirectly because the 
chaotic dynamical interactions between protostars and ejections of 
stars and brown dwarfs from the cloud produce different gravita- 
tional potentials and move the gas distributions in different ways. 

BBB2003 also presented snapshots of the evolution in two 
lower-mass dense cores that produced 7 and 5 objects, respectively. 
In each core, all but one of these objects formed via the fragmenta- 
tion of a disc surrounding an protostar. With radiation hydrodynam- 
ics, these two cores produce 4 and 1 objects, respectively. These 
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Figure 5. Time of formation and mass of each star and brown dwarf at the 
end of original BBB2003 calculation (top) and the radiation hydrodynami- 
cal calculations BBB2003 RT5 (centre) and BBB2003 RT0.5 (bottom). Ra- 
diative feedback dramatically decreases the number of objects formed, pri- 
marily by suppressing the continued star formation in the main dense core 
that occurs in the original calculation. Objects that form in the main dense 
core are denoted by blue lines. Objects that form in the lower-mass sec- 
ond and third dense cores are denoted by green and red lines, respectively. 
Objects that are still accreting significantly at the end of the calculation 
are represented with vertical arrows. The horizontal dashed line marks the 
star/brown dwarf boundary. Time is measured in terms of the free-fall time 
of the initial cloud (top) or years (bottom). 



two dense cores are shown at the end of the calculation in Figures 
[3]for comparison with BBB2003. 

Although radiative feedback strongly reduces the fragmenta- 
tion on small length scales, this does not prohibit the formation of 
binary and multiple systems. With only 13 objects produced by the 
calculation (Table[T}, it is impossible to discuss multiplicity statis- 
tics. However, we note that when the calculation is stopped there is 
one system of six objects (two 2-AU binaries orbiting each other at 
20 AU and this system is orbited by a 27-AU binary at 430 AU), 
one quintuple system (a 17-AU binary with companions orbiting 
at 65, 234, and 13000 AU), and two single stars. As was found by 
|Bate e t al. (2002), although two objects can form well separated 
from each other, the combination of dynamical interactions, gas 
accretion, and interactions with circumbinary and circum-multiple 
discs is very effective at producing multiple systems and even close 
binaries. For example, the second and third objects to form in the 
calculation are initially separated by more than 1600 AU and yet 
at the end of the calculation they comprise one of the two 2-AU 
binaries. 

In Figure [5] we plot the final mass of each object versus 
the time of its formation (i.e. the time that a sink particle was 
inserted) for the original BBB2003 calculation (top panel), the 
BBB2003 RT5 calculation (centre panel) and the BBB2003 RT0.5 
calculation (lower panel). The second main difference between the 
barotropic calculation and the radiation hydrodynamical calcula- 
tions is that after the first burst of star formation in the main dense 
core (t — 1.03 — 1.13 tg for each of the simulations, which is 
approximately the dynamical timescale of the dense core) there 
is a second burst of star formation in the barotropic calculation 
(t = 1.25 — 1.35 tff) but not in the radiation hydrodynamical sim- 
ulations (BBB2003 RT5 has a weak burst at t « 1.27 tg , but only 
one further object is formed in the main dense core of BBB2003 
RT0.5). Radiative feedback is responsible for almost shutting off 
the production of new stars in the main dense core. Whereas stel- 
lar feedback is often considered one of the main agents in termi- 
nating the star formation process, the feedback invoked is usually 
that from massive stars, not low-mass stars (e.g. the formation of 
HII regions, strong winds, and supernova explosions). The calcu- 
lations presented here show that even the radiative feedback from 
low-mass star formation is enough to almost terminate the produc- 
tion of new stars in dense molecular cores with masses ^ 10 Mq. 
It does not destroy the dense molecular core or stop the existing 
protostars from accreting more gas, but it essentially turns off the 
production of new objects by raising the temperature of the bulk of 
the dense core so that few new gravitationally-unstable condensa- 
tions are formed. 



3.1.2 Sink particles with r acc = 5 AU 

In Figure|6] we present snapshots of the evolution of the main dense 
core from the BBB2003 RT5 calculation (i.e. the same as above, but 
with sink particle accretion radii 10 times larger). These snapshots 
should be compared with those in Figure [2] The initial evolution is 
very similar to that obtained with r acc = 0.5 AU but it can be seen 
that the gas temperatures are slightly lower due to the fact that less 
accretion luminosity is injected into the calculation. From the snap- 
shots and Figure|5]it can be seen that up until t ~ 1.12 tg the two 
calculations evolve in a similar manner. However, the 5th object to 
form in BBB2003 RT5 is prevented from forming in the BBB2003 
RT0.5 calculation due to the slightly higher gas temperatures in the 
central region visible in the t = 1.08 and 1.12 tg panels of Fig- 
ures|2]and|6] In BBB2003 RT5, this object is formed only 24 AU 
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Figure 6. The star formation in the main dense core of the BBB2003 RT5 calculation. These images are identical to those from the BBB2003 RT0.5 calculation 
in Figure[2]except that the calculation here uses sink particle accretion radii that are an order of magnitude larger (5 AU rather than 0.5 AU). The early evolution 
of the cloud (top panels) is relatively independent of the sink particle parameters, but it can be seen that the gas near the prototars is slightly cooler with the 
larger accretion radii because a larger fraction of the accretion luminosity is neglected. Later (lower panels), the two calculations slowly diverge. There is 
slightly more fragmentation in the larger accretion radius calculation (but still substantially less than in the original barotropic BBB2003 calculation). 



from an existing protostar. Subsequently in this central region, only 
one of the 8th, 9th, and 10th objects to form in BBB2003 RT5 (all 
within 350 AU of existing protostars) manages to collapse in the 
BBB2003 RT0.5 calculation. Thus, it can be seen that neglecting 
the accretion luminosity emitted from within 5 AU of the protostars 
does increase the amount of fragmentation (although even this lim- 
ited radiative feedback still reduces the number of objects formed 
by more than a factor of 3 compared to the barotropic equation 
of state). Whether or not the 0.5 AU accretion radii calculation is 
converged in the sense that using smaller accretion radii would not 
change the number of fragments is, of course, hard to tell. However, 
the small accretion radii calculation only forms two fewer objects 
than the large accretion radius calculation with the enhanced radia- 



tive feedback (Table [TJ, and only two protostars form closer than 
100 AU from existing protostars in the small accretion radii calcu- 
lation (at 23 and 30 AU separations). Thus, modelling the accretion 
flows right down to the surfaces of the protostars is unlikely to de- 
crease the number of objects much further. 

Beyond t = 1.14 tg, the calculations diverge because of the 
different numbers of objects and their chaotic dynamical interac- 
tions and ejections. Interestingly, this even affects the star forma- 
tion in the two lower-mass dense cores ~ 2 x 10 4 AU from the main 
dense core. In the smaller accretion radii calculations each of these 
two dense cores begin forming their protostars w 10, 000 years 
earlier than in the larger accretion radii calculations. This is pre- 
sumably due to slight perturbations of the lower-mass dense cores 
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Figure 7. Histograms giving the differential initial mass functions of the 13, 15, and 17 stars and brown dwarfs that had been produced by the end of the 
BBB2003 RT0.5, BBB2003 RT5, and BB2005 RT0.5 calculations, respectively. The single hashed region gives objects that have stopped accreting while the 
double hashed region gives those objects that are still accreting. Parameterisations of the observed IMF by Salpeter 1 1955 1, Kroupa 1 2001 1 and Chabrier 1 2003 1 
are given by the magenta line, red broken power law, and black curve, respectively. 
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Figure 8. The cumulative IMFs produced by all 5 calculations discussed 
in this paper. The previously published IMFs from BBB2003 and BB2005 
using a barotropic equation of state are given by the thin solid line and thin 
dashed line, respectively. The radiation hydrodynamical calculations pre- 
sented here are BBB2003 RT0.5 (thick solid line), BB2005 RT0.5 (thick 
dashed line) and the large accretion radius BBB2003 RT5 calculation (thin 
dot-dashed line). The vertical long-dashed line denotes the boundary be- 
tween brown dwarfs and stars. It is clear that the radiation hydrodynamical 
calculations produce IMFs with a larger characteristic mass and far fewer 
brown dwarfs and low-mass stars than the original barotropic calculations. 
Furthermore, whereas BBB2003 and BB2005 showed a clear dependence 
of the characteristic stellar mass on the initial Jeans mass of the molecular 
clouds (BB2005 began with a denser cloud with a Jeans mass 3 times lower 
that produced a median stellar mass 3.04 times lower than BBB2003), when 
radiative feedback is included there is no significant dependence of the IMF 
on cloud density and the initial Jeans mass. A Kolmogorov-Smirnov (K-S) 
test on the BBB2003 RT0.5 and BB2005 RT0.5 distributions gives a 99.97% 
probability that the two IMFs were drawn from the same underlying distri- 
bution (i.e. they are statistically indistinguishable). By comparison, the two 
IMFs from the original barotropic calculations had only a 1 .9% probability 
of being drawn from the same underlying distribution. 



caused by the differing gravitational potentials due to the different 
evolutions of the main dense core. 



3.1.3 The initial mass functions 

The effect of radiative feedback in terminating the production of 
new objects within the dense cores after one dynamical timescale 
and inhibiting the fragmentation of discs and filaments near ex- 
isting protostars has a tremendous effect on the number of ob- 
jects formed and the final distribution of stellar masses. Table [I] 
summaries the numbers of stars and brown dwarfs formed, their 
combined mass, and their mean and median masses. The origi- 
nal barotropic BBB2003 calculation produced 50 stars and brown 
dwarfs in 1.40 t a . However, BBB2003 RT0.5 only produced 13 ob- 
jects in the same time and even BBB2003 RT5 with less accretion 
luminosity from the protostars produced only 15 objects. Thus, the 
inclusion of radiative feedback has cut the number of objects pro- 
duced by a factor of « 4. In addition, whereas the original calcula- 
tion produced a similar number of stars and brown dwarfs the ratio 
of brown dwarfs to stars is 1:3 for BBB2003 RT5. For BBB2003 
RT0.5, the ratio is less than 1:5 and both objects with brown dwarf 
masses are still accreting when the calculation is stopped. 

The much lower fraction of brown dwarfs is due to both the 
inhibiting of the fragmentation of discs and nearby filaments (be- 
cause such objects are frequently ejected through dynamical inter- 
actions, terminating their accretion before they have been able to 
accrete much mass) and the suppression of new objects formed 
in the dense cores after a dynamical time. In the latter case, it 
can be seen in Figure [5] that there is a higher fraction of brown 
dwarfs amongst objects formed later in the barotropic calculation 
than those formed earlier (top panel). Objects that form later must 
compete with the higher-mass protostars for the available gas. Usu- 
ally they lose, either being dynamically ejected from the dense core 
or at least having their velocities increased relative to the gas so that 
their accretion rates drop (Bondi-Hoyle accretion is proportional to 
l/v 3 ). Ejections producing brown dwarfs and low-mass stars still 
occur, but they are much less common with the inclusion of radia- 
tive feedback than they were in the barotropic calculation. 

Although the number of objects is decreased by the inclusion 
of radiative feedback, the amount of gas that has been converted 
into stars at t = 1.40 ts is actually about 15% greater compared to 
the barotropic calculation (see Table[T](. 

The overall result of all of these effects is that the characteris- 
tic mass of the IMF moves to higher masses with the inclusion of 
radiative feedback and fewer brown dwarfs and low-mass stars are 
produced. Comparing BBB2003 RT0.5 with BBB2003, the mean 
and median masses have increased by a factor of 4.4 to ~ 0.5 Mq 
and ~ 0.3 Mq, respectively (Table [TJ. 
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Unfortunately, with so few objects it is not really possible to 
plot meaningful differential IMFs, but for completeness we include 
these in Figure[7]in a similar manner to the original BBB2003 pa- 
per. The IMFs are compared with the parameterisations of the ob- 
served IMF given by |Chabrier|p003) , |Kroupa| ( |200T) , and |Salpeter| 
( [1955} . 

In Figure [8] we plot the cumulative IMFs produced by the 
original barotropic calculations and the three radiation hydrody- 
namical simulations presented here. The BBB2003-type calcula- 
tions are given using the solid lines (thin for the barotropic equation 
of state and thick for BBB2003 RT0.5) and the dot-dashed line for 
the larger accretion radius BBB2003 RT5 calculation. The shift of 
the IMF to higher masses with the inclusion of radiative feedback 
is clear. A Kolmogorov-Smirnov (K-S) test comparing the original 
calculation with BBB2003 RT0.5 gives only a 1.7% probability that 
the two IMFs are drawn from the same underlying population (i.e. 
they are different). We note that a K-S test comparing BBB2003 
RT5 with BBB2003 RT0.5 shows they are indistinguishable (there 
is a 35% probability they are drawn from the same population). 

3.2 BB2005 initial conditions 

We turn now to the second set of initial conditions, those from 
BB2005 which are identical to those of BBB2003 except that the 
initial cloud has 9 times the density due to a smaller radius, and a 
higher Mach number so that for both types of initial conditions the 
kinetic and gravitational potential energies have the same magni- 
tude (see Table[TJ. Since the cloud has the same initial temperature, 
the mean thermal Jeans mass is a factor of 3 lower (1/3 M© rather 
than 1 Mo). In the original barotropic calculations, the lower initial 
Jeans mass resulted in a proportional lowering of the median stellar 
mass of the objects produced, showing that the characteristic stellar 
mass is set by the initial Jeans mass in such simulations. 

The radiation hydrodynamical version of the BB2005 calcu- 
lation (BB2005 RT0.5) uses sink particles with the more accurate 
small accretion radii (r acc = 0.5 AU). Again, we do not begin the 
radiation hydrodynamical simulation from the initial conditions be- 
cause the early evolution of the cloud is essentially isothermal. We 
begin calculation BB2005 RT0.5 from the last dump of BB2005 
before the density exceeded 9 x 1CP 16 g cm~ 3 at t = 0.692 tg. In 
the original BB2005 calculation the first sink particle was inserted 
at t — 0.824 iff, some 8.4 x 10 3 years later. Using radiation hy- 
drodynamics again delays the insertion of the first sink particle (for 
the same reason as given in Section [3~T) until t = 0.850 ta. 

In the righthand panels of Figure ^ we present snapshots 
of the global evolution of the BB2005 RT0.5 calculation. Fig- 
ure [9] shows the evolution of the main dense core from the 
BB2005 RT0.5 calculation in much greater detail than Figure [TJ 
Again, these are shown at the same times as the equivalent fig- 
ures in the original BB2005 paper to allow direct comparison. 
Again, an animation comparing the original calculation with the 
radiation hydrodynamical calculation can be downloaded from 
http://www.astro.ex.ac.uk/people/mbate~ 

As with the BBB2003-type initial conditions, the barotropic 
and radiation hydrodynamical calculations diverge quickly on 
small scales. In the barotropic calculation, the first two objects form 
from two dense clumps separated by approximately 100 AU, one 
object forming 700 years after the other. This binary accretes a 
massive circumbinary disc which later fragments to form a third 
object. In the radiation hydrodynamical simulation, the accretion 
luminosity from the first object inhibits the formation of the second 
object and the massive disc which subsequently forms around this 



single protostar is too hot to fragment. Thus, whereas three objects 
were formed in this small region of the original simulation, the ra- 
diation hydrodynamical simulation forms only a single protostar. 
This pattern continues throughout the early evolution - fragmen- 
tation on length scales J> 300 AU is unaffected by the accretion 
luminosity emitted by existing protostars, but most of the fragmen- 
tation on smaller scales (of either nearby dense filaments or massive 
circumstellar discs) is stopped. Figure]?] shows that more than 3/4 
of the objects form within 300 AU of an existing protostar in the 
original barotropic BB2005 calculation (thin dashed line) whereas 
this fraction is reduced to less than 1/3 with radiative feedback 
(thick dashed line). As with the BBB2003-type initial conditions, 
the radiative feedback eventually begins to alter structures on larger 
scales indirectly because of the differing numbers of protostars and 
their dynamical interactions (i ^ 1.20 iff). 

As discussed above for the first type of initial conditions, the 
reduction of small-scale fragmentation due to the radiative heat- 
ing does not stop the production of binary and multiple objects. At 
the end of the calculation there are 2 single very-low-mass objects 
(masses < 0.1 Mq), a triple system consisting of a 150- AU binary 
with a wide (8400-AU) companion, and a small bound cluster of 12 
objects, including three binaries with semi-major axes of less than 
3 AU, one of which has a close third companion with a semi-major 
axis of 20 AU. 

In Figure [10] we plot the final mass of each object versus the 
time of its formation (i.e. the time that a sink particle was inserted) 
for the original BB2005 calculation (top panel) and the BB2005 
RT0.5 calculation (lower panel). As with the first type of initial con- 
ditions, the radiative feedback terminates the production of new ob- 
jects in the main dense core well before the simulation is stopped. 
By contrast, in the barotropic calculation, production of stars and 
brown dwarfs in the main dense core continues until the calcula- 
tion is stopped, with only a brief pause at t — 1.22 — 1.30 tg. With 
radiative feedback the main dense core ceases production of new 
objects at t = 1.20 because the gas in the main dense core has 
been heated by the embedded protostars. 

In the original BB2005 calculation, star formation proceeded 
in three lower-mass dense cores in addition to the main dense core. 
This is also true of the radiation hydrodynamical simulation. Two 
of these dense cores (with masses 0.2 — 0.3 Mq when the sink par- 
ticles are inserted) produce single protostars in both the barotropic 
and radiation hydrodynamical calculations. The remaining dense 
core (with a mass of « 1 M Q when star formation begins) formed 
12 objects in the barotropic calculation due to a combination of disc 
fragmentation and filament fragmentation. With radiation hydrody- 
namics, this core only produces 2 objects due to fragmentation of a 
filament. 



3.2.1 The initial mass function 

As with the lower-density cloud, the radiative feedback dramati- 
cally decreases the number of objects formed (see Table [JJl. The 
original BB2005 calculation produced 79 stars and brown dwarfs 
in 1.40 iff, whereas with radiation hydrodynamics only 17 objects 
form. Whereas the barotropic calculation produced three times as 
many brown dwarfs as stars, this ratio is more than inverted when 
radiative feedback is included, with more than 4 times as many stars 
as brown dwarfs (including one brown dwarf that has had its ac- 
cretion terminated). Again, this lower proportion of brown dwarfs 
is due to the decrease in disc fragmentation and the collapse of 
nearby filaments, as well as the cessation of fragmentation in the 
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Figure 9. The star formation in the main dense core of the BB2005 RT0.5 calculation. The first object forms at t = 0.850 tg. Large gaseous filaments 
collapse to form single objects and multiple systems. These objects fall together to form a small group. Dynamical interactions within the group eject a few 
objects. Radiative feedback from the accreting protostars heats the gas in the dense core. Each panel is 0.025 pc (5150 au) across. Time is given in units of 
the initial free-fall time of 6.34 X 10 4 yr. The red-yellow-white panels show the logarithm of column density, N, through the cloud, with the scale covering 
— 0.2 < log N < 2.5 with N measured in g cm" 2 . The blue-red-yellow-white panels show the logarithm of mass weighted temperature, T, through the cloud 
with the scale covering 9 — 300 K. 



main dense core at t — 1.20tg due to the heating of the core by the 
embedded protostars (Figure[T0}. 

The amount of gas converted into stars and brown dwarfs at 
the end of the calculations is similar in the barotropic and radiation 
hydrodynamical calculations (only differing by 4%: Table [TJ, but 
as with the BBB2003-type initial conditions the large decrease in 
the number of objects produced moves the characteristic mass of 



the IMF to much higher masses. The mean and median masses are 
increased by factors of 4.5 and 15, to « 0.5 Mq and « 0.3 Mq, 
respectively. 

The differential and cumulative IMFs are displayed in Figures 
[7] (right panel) and [5] (thick dashed line), respectively, in a similar 
manner to those given in the original BB2005 paper. The barotropic 
cumulative IMF is also plotted in Figure [8] (thin dashed line), 
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Figure 10. Time of formation and mass of each star and brown dwarf at the 
end of original BB2005 calculation (top) and the radiation hydrodynami- 
cal calculation BB2005 RT0.5 (bottom). Again, radiative feedback does not 
alter the initial phase of star formation too much, but most of the contin- 
ued star formation that occurred in the original calculation is inhibited by 
the radiative feedback. Objects that form in the main dense core are de- 
noted by blue lines. Objects that form in the lower-mass second, third, and 
fourth dense cores are denoted by green, red, and magenta lines, respec- 
tively. Objects that are still accreting significantly at the end of the calcula- 
tion are represented with vertical arrows. The horizontal dashed line marks 
the star/brown dwarf boundary. Time is measured in terms of the free-fall 
time of the initial cloud (top) or years (bottom). 



clearly displaying the increase in the typical stellar mass. A K-S test 
comparing the original barotropic calculation with BB2005 RT0.5 
gives only an 0.0004% (i.e., 1/250,000) probability that the two 
IMFs are drawn from the same underlying population (i.e. they are 
distinctly different). 



4 DISCUSSION 

4.1 The ratio of stars to brown dwarfs 

The observed initial mass function within the local region of our 
Galaxy is now reasonably well constrained, at least down to « 0.03 
Mq. Surprisingly, there is little evidence for variation of the IMF 
amongst star-forming regions, open clusters, and even globular 
clusters, leading to the idea that the IMF may be universal (at least 



for metallicities greater than ~ 1/1000 solar). Two frequently used 
parameterisations of the observed IMF are given by Kroupa ( 2001 ) 
and Chabrier (2003) (see Figure [7](. In particular, it is now gener- 
ally accepted that the number of stars is larger than the number of 
brown dwarfs jChabrier|2003l |Greissl et al.|2007| |Luhman|2007l 
|Andersen et al.|20 08). Andersen et al. (2008) find that the ratio of 
stars with masses 0.08 — 1.0 Mq to brown dwarfs with masses 
0.03 - 0.08 M is N(0.08 - 1.0)/JV(0.03 - 0.08) » 5 ± 2. 

Previous hydrodynamical simulations using a barotropic equa- 
tion of state have consistently produced more brown dwarfs than 
stars (BBB2003, BB2005, |Bai^|2005l |2009"1 >. In the earlier cal- 
culations, this might have been explained away by small number 
statistics (the earlier calculations each produced fewer than 100 ob- 
jects). However, Bate (20091 recently performed a barotropic cal- 
culation of the collapse and fragmentation of a 500 Mq cloud that 
produced 1254 stars and brown dwarfs. While some of the stellar 
properties (e.g. the multiplicity as a function of primary mass) were 
in good agreement with observations, this calculation produced at 
least 50% more brown dwarfs than stars. The large number of ob- 
jects produced leaves no doubt as to the statistical significance of 
this result. 

For this reason, the decrease in the fraction of brown dwarfs 
to stars provided by the inclusion of radiative feedback is welcome. 
Averaging over the two small accretion radii calculations presented 
here, the ratio of stars to brown dwarfs is found to be <^ 5. Although 
the statistics are poor, this is at least in reasonable agreement with 
observations. Obviously the next step is to perform simulations of 
more massive molecular clouds in order to improve the statistics as 
was done by Bate (2009) for the barotropic calculations. 

Magnetic fields may also affect the mass function. Recently, 
|Price & Bate ( 2007 1 showed that stronger magnetic fields generally 
inhibit disc formation and binary formation (see also Hennebelle & 
|Fromang|2008"l|Hennebelle & Teyssier|2008| l. |Pri"ce & BateH20"08l > 
ran star cluster formation simulations similar to BBB2003, but with 
magnetic fields. They found that the extra pressure support pro- 
vided by magnetic fields generally decreased the rate of star forma- 
tion and the importance of dynamical interactions between objects. 
This leads to stronger magnetic fields producing a decrease in the 
ratio of brown dwarfs to stars (though again the total numbers of 
objects formed in the calculations were small, ranging from 15 to 
69). Finally, |Qffner, Klein & McKee|f2008fr report that simulations 
of driven turbulence produce fewer low-mass objects than simula- 
tions with decaying turbulence. 



4.2 Comparison with the Taurus-Auriga star-forming region 

The Taurus-Auriga star-forming region is a large but diffuse molec- 
ular cloud spanning 30 pc and containing 2 x 10 4 Mq of molec- 
ular gas (Golds mith et al.| |2008). Although the entire complex is 
much larger and more massive than the simulations presented here, 
many of the recently formed stars in Taurus are contained within 
6 or 7 groups each containing ~ 10 members and measuring ~ 2 
pc across ( jGomez et al.|1993] >. These groups may have originated 
from dense cores similar to the two main dense cores modelled in 
the calculations presented here which form 8 and 13 objects. Once 
dispersed, these groups would appear very similar to the Taurus 
groups. The IMF of the Taurus cloud (Lu hman et al.|2 003 1 is also 
similar to the IMFs presented in Figure [7] in the sense that they 
both seem to contain an excess of ~ 1 Mq stars when compared to 
parameterisations of a universal IMF. 

Although the Taurus IMF and the simulation IMFs are con- 
sistent with a universal IMF due to the small number statistics, it 
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may also be that low-density distributed star-forming regions that 
form stars in small groups rather than large clusters do preferen- 
tially form stars with a masses ~ 1 Mq. In the simulations pre- 
sented here, the small groups form in dense molecular cores with 
masses of « 10 Mq that fragment into « 10 objects. With a few 
objects dynamically ejected as brown dwarfs or low-mass stars be- 
fore they have accreted much of the available mass, the remaining 
objects can each accrete to m 1 Mq before the gas reservoir is 
exhausted. This leaves a stellar mass distribution that is biased in 
favour of solar-type stars, with a few lower mass members. Such 
an explanation for the IMF in Taurus was first presented by |Good-] 
win, Whitworth & Ward-Thompson ( 2004 1 based on hydrodynam- 
ical simulations of isolated dense cores. If this model is correct, the 
simulations presented here may be best compared with the small 
Taurus stellar groups. 



|Nayakshin|(2006l > proposes that radiative feedback from low-mass 
protostars forming in a gravitationally-unstable disc surrounding a 
supermassive black hole may heat the disc enough to stop further 
fragmentation. Trapped in the disc, the few protostars that did man- 
age to form before the disc was stabilised by the protostellar heat- 
ing would then accrete the remaining gas becoming very massive. 
Nayakshin proposes that this may explain the apparent top-heavy 
mass function of the stars orbiting the Sgr A* in the centre of our 
Galaxy (Nayakshin & Sunyaev 20051. Although this case is very 
different to local star formation, it is a case where radiative feed- 
back, even from low-mass protostars, may play a substantial role 
in setting stellar masses. In what follows, we show that radiative 
feedback may not just set the characteristic stellar mass in this ex- 
otic case - it may also be responsible for setting the characteristic 
stellar mass in typical star formation. 



4.3 The characteristic mass of the IMF 

The calculations of BBB2003, BB2005, and [Bate] ( |2003) tested the 
dependence of the IMF obtained from barotropic hydrodynamical 
simulations on various changes to the initial conditions and metal- 
licity of molecular clouds. The conclusion from these papers was 
that the characteristic mass of the IMF is set by the initial Jeans 
mass of the molecular clouds. Other similar hydrodynamical cal- 
culations, that did not resolve brown dwarfs or close binaries, have 
led to similar conclusions (Klessen et al. 1998; Klessen & Burkert 
|2000| |Jappsen et al.|2005| |Bonnell et al.||2006| l. Current thinking 
is that this initial Jeans mass may be set by the thermodynamics 
of molecular gas at the transition from atomic line cooling to dust 
cooling which sets a characteristic density and Jeans mass ( |Larson| 
|1985||2005| . Again, this appears to be backed up by hydrodynam- 
ical simulations showing that the peak of the IMF scales linearly 
with a change in the density at which this transition occurs jjappsen| 
et al. 2005 , Bonnell et al. 20061. Even other models of the origin 
of the IMF predict that the characteristic mass of the IMF should 
scale linearly with the typical Jeans mass. For example, the theory 
of |Padoan & Nordlund| < |2002"l > which proposes that the IMF origi- 
nates from the mass spectrum of dense cores formed by supersonic 
turbulence predicts that the characteristic mass of the IMF scales 
linearly with the Jeans mass of the molecular cloud but also de- 
pends on the Mach number of the turbulence. 

Recently, Elmegreen et al. ( 2008 1 proposed that the reason the 
IMF appears to be a universal function in the local Universe is be- 
cause this characteristic thermodynamic Jeans mass in molecular 
clouds is relatively insensitive to initial conditions such as molec- 
ular cloud density, the local radiation field, etc. However, for the 
dense cores where stellar groups form this theory implies that tem- 
perature at dust-gas coupling increases with density as n 1//2 . If any- 
thing, observations tend to show the temperature in dense molec- 
ular gas decreases with increasing density. This theory also means 
that the IMF is universal because of the thermodynamics of gas and 
dust, not due to the star formation process itself. It would be more 
elegant if star formation somehow regulated itself to provide a near 
universal IMF in the local Universe. 

Given the wide range of possible theories for the origin of the 
IMF (see Section[TJ, little attention has been paid to the role of ra- 
diative feedback. Dynamical feedback processes such as winds and 
outflows have been considered ( |Shu et al.|1988||Silk|1995||Adarrisl 
& Fatuzzo 1996), but the effect of radiative feedback on stellar 
masses is usually only considered for massive stars in the context 
of setting a maximum stellar mass. One exception is star formation 
in massive accretion discs surrounding supermassive black holes. 



4.3.1 The effect of radiative feedback on the characteristic mass 

In the barotropic calculations of BBB2003 and BB2005, the density 
of the initial molecular cloud was 9 times higher in the latter cal- 
culation leading to an initial mean thermal Jeans mass three times 
lower in the denser cloud. The mass distributions of the stars and 
brown dwarfs produced by these calculations had median masses 
that scaled almost exactly with this change in the initial Jeans mass 
- the median mass was 3.04 times smaller for the stellar cluster 
produced by the denser cloud. K-S tests performed on the two dis- 
tributions showed that there was less than a 1.9% chance that they 
were drawn from the same underlying mass function. 

However, as seen above, when the calculations are repeated 
with radiative feedback the resulting IMFs have a significantly 
larger characteristic mass and, more importantly, they are indistin- 
guishable from each other (Figure [Hf despite the different initial 
mean thermal Jeans masses of the clouds. A K-S test comparing 
the IMFs of the BBB2003 RT0.5 and BB2005 RT0.5 calculations 
gives a 99.97% probability they are drawn from the same popula- 
tion. This implies that radiative feedback during the star formation 
process substantially weakens the dependence of the IMF on the 
initial Jeans mass in the molecular cloud and may be responsible 
for producing a nearly universal IMF. 

The general idea of how radiative feedback may increase the 
characteristic mass of the IMF is relatively simple to understand. 
The Jeans length scales with temperature and density in a molecular 
cloud as 

Aj oc r 1/2 p - 1/2 , (i) 

so the Jeans mass scales as 

Afj oc pAj oc T 3/ V 1/2 . (2) 

Thus, for the two types of initial conditions here which have the 
same temperature but a density that differs by a factor of 9, the 
initial Jeans mass and Jeans length are 3 times smaller in the denser 
cloud. When the gas collapses, the fragments (protostars) that form 
are roughly separated by the Jeans length so that in the denser cloud 
they are closer together and, therefore, each has a smaller reservoir 
from which to accrete. Essentially, a molecular cloud of mass M c 
is broken into M c /Mj objects that, on average, accrete the mass 
available in the reservoir contained within a Jeans length and, thus, 
their characteristic mass is ~ Mj (some objects accrete more and 
some less due to competitive accretion, but the characteristic mass 
is « Afj). 

However, radiative feedback changes this situation. When the 
first fragment forms, this protostar heats the gas around it. Nearby 
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gas, which would otherwise have collapsed soon after the first ob- 
ject to form additional objects, may now be prohibited from col- 
lapsing because it is hotter and that gas is instead accreted by the 
already existing protostar. Thus, the distance between neighbouring 
protostars is increased because of the radiative feedback and those 
that do form each have a larger reservoir of gas to accrete. In other 
words, the effective Jeans length and Jeans mass (the mass that will 
end up in that object rather than collapsing to form a neighbour- 
ing object) has increased. This explains the general increase in the 
characteristic mass. 



4.3.2 Dependence of the characteristic mass on initial conditions 

Understanding how the characteristic mass of the IMF may depend 
on variations in the initial conditions (e.g. the cloud's mean density 
or temperature) when radiative feedback is included is somewhat 
more complicated. 

Let us begin by assuming that the mean density of the molec- 
ular cloud is low enough that the cloud is optically thin to in- 
frared radiation. At a typical temperature of 10 K for solar metal- 
licity molecular gas the mean dust absorption opacity is k w 0.01 
cm 2 g" 1 so that for molecular hydrogen densities of 10 6 cm -3 an 
optical depth of unity is reached after w 10 pc, much larger than the 
typical distance between protostars in a typical star-forming region. 
Even taking a higher density of 10 7 cm' 3 and a higher temperature 
of 30 K with a corresponding mean opacity of k — 0.1 cm 2 g _1 
this length scale is 0.1 pc. For simplicity, and because the radiation 
hydrodynamical calculations presented in this paper have been per- 
formed using grey radiative transfer, we also assume that the dust 
absorption opacity is independent of wavelength. The more general 
case, in which the wavelength dependence of the dust absorption 
opacity scales as g a b s oc A - * 3 , is relatively easy to derive but the 
extra complication distracts from the main point of the discussion 
below. 

Under these assumptions, the temperature, T, of the gas and 
dust at radius, r, from a (spherically symmetric) protostar of lumi- 
nosity, L«, is given by 

L* = 47ro-r 2 T 4 , (3) 

where a is the Stephan-Boltzman constant. 

Now, a convenient definition of the Jeans length, A.j, is that it 
is the radius at which the sound speed of the gas equals the escape 
velocity of the mass enclosed within this radius 

GM 



— T 
/' 



(4) 



where c s is the sound speed, TZ is the gas constant, /i is the mean 
molecular weight, and the Jeans mass is 
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On scales smaller than Aj, the gas is supported by pressure against 
collapse, while on larger scales the gas is unstable to collapse. For 
the usual definition of the Jeans length and mass, the sound speed 
is taken to be a constant and by combining equations [4] and [5] we 
obtain equations[T]and|2]. 

However, for a gas cloud that is internally heated by an em- 
bedded protostar we can use equations|3]|4] and|5]to obtain 
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where we have used A c fi and M e e to differentiate the effective 
Jeans length and mass due to radiative feedback from the stan- 
dard (isothermal) definitions, Aj and Mj. Thus, the effective Jeans 
mass in a molecular cloud containing accreting protostars is signif- 
icantly less dependent on the density of the cloud than the usual 
Jeans mass, Mj (which scales oc p" 1 / 2 ). This dependence of the 
effective Jeans mass on density is so weak that it is consistent with 
the results of the radiation hydrodynamical simulations presented 
here in which we find no significant variation of the characteris- 
tic mass of the IMF with density. Simulations that formed a much 
greater number of objects would be needed to detect such a weak 
dependence. We also note that the dependence on density in equa- 
tion [7] is even weaker than that given by the theory of |Elmegreen| 
|et aL ^2008} , and it has the opposite sign (i.e. the characteristic mass 
decreases with increasing density whereas Elmegree n et al.| find 
that it should increase with increasing density). Finally, as men- 
tioned above, we have neglected the wavelength dependence of the 
dust absorption opacity. Including this alters the powers of density 
and luminosity that appear in equations|6]and|7]slightly. For exam- 
ple, the scaling of the effective Jeans mass with density becomes 
p-!/(5-/3) so tnat f or with (3—1 the effective Jeans mass 
scales as p~ 1//4 . This is still substantially weaker than the usual 
scaling of Jeans mass with cloud density. 

Evaluating this equation is not straightforward because of the 
dependence on the protostellar luminosity which depends on the 
mass, radius and accretion rate of the protostar. However, because 
the dependence on the protostellar luminosity is quite weak, we can 
put in estimated numbers to check that it gives a realistic character- 
istic mass. The luminosity of an accreting protostar is 



GM,M* 



(8) 



where M* and i?» are the mass and radius of the protostar, respec- 
tively. In the radiation hydrodynamics calculations presented here, 
the mean accretion rate of the objects is approximately the same for 
both calculations: 1.5 x 10~ 5 M n yr" 1 for BBB2003 RT0.5 and 

2.5 x 10~ 5 Mr. 



1 for BB2005 RT0.5 (Figure 
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These rates 



are calculated as the mean of the time-averaged accretion rates of 
all the objects in a calculation, where the time-averaged accretion 
rate is the mass of the object at the end of the calculation divided 
by the time between its formation and the end of its accretion or 
the end of the calculation, which ever comes first. Assuming the ra- 
diative feedback typically comes from a 0.1 Mq protostar of radius 
2 R.0 accreting at a rate of 1 x 10 M© yr -1 gives a luminosity 
of 150 Lq. The effective Jeans mass for the BBB2003-type initial 
conditions with an initial cloud density of 1.2 x 10" 19 g cm" 3 can 
then be written 



0.5 



1.2 x 10" 



5 / L* 



V 150 L e 



M . (9) 



This is a reasonable value for the characteristic mass of the IMF. 

To try and proceed beyond this rough estimate, we need to 
make further simplifying assumptions about the luminosity of the 
protostars. The protostellar accretion rates may be expected to scale 
as ~ c 3 /G. However, even for an isothermal, non-singular spheri- 
cal cloud the accretion rate usually begins much larger than c 3 /G 
and declines with time (Foster & Chevalier 1993). Here the situa- 
tion is even more complex because the sound speed is that of the 
internally heated cloud which is a function of the radius from the 
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Figure 11. The time-averaged accretion rates of the objects versus their fi- 
nal masses from the BBB2003 RT0.5 (top) and BB2005 RT0.5 (bottom) 
calculations. The accretion rates are calculated as the final mass of an ob- 
ject divided by the time between its formation and the termination of its 
accretion or the end of the calculation. The mean accretion rates for the 
two calculations are 1.5 X 10~ 5 Mq yr" 1 and 2.5 X 10~ 5 Mq yr _1 , re- 
spectively. These values are statistically indistinguishable from each other 
given the small numbers of objects and the fact the the dispersions of the 
accretion rates are larger than the difference between the two calculations. 
The horizontal solid lines give the mean accretion rates. Objects still accret- 
ing when the calculations were stopped are denoted with horizontal arrows. 
The accretion rates are given in M-@/tg on the left-hand axes and Mo/yr 
on the right-hand axes. The vertical dashed line marks the star/brown dwarf 
boundary. 



protostar rather than simply the initial temperature of the cloud. On 
the other hand, the timescale for the accretion of the envelope must 
still be ~ A c ff / c s where the relevant sound speed is that at distance 
A c ff from the protostar. Therefore, we can write 
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Again using equation[3]to eliminate T and then using equation[8]to 
solve for L, we obtain 
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Inserting this into equation[6]we get 
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which finally gives 
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mass is that the dependence on the initial density of the cloud has 
been almost eliminated and all that remains is a dependence on the 
protostellar mass to radius ratio. This result should be treated with 
caution, simply because of the number of simplifying assumptions 
that have been involved. Furthermore, it is quite likely that the mass 
to radius ratio of the protostar may itself depend on the mass accre- 
tion rate in which case equation [13] still contains a dependence on 
the initial density of the molecular cloud via equation[l0] However, 
the main point to take away from this discussion is that in both 
equations [7] and [T3] the effective Jeans mass is at most only weakly 
dependent on the initial density of the molecular cloud which is 
consistent with our numerical results. 

Finally, we note that we have assumed in the above discussion 
that the protostar's luminosity provides the dominant contribution 
to the temperature at distances Si A e fi from the protostar, as op- 
posed to the background temperature. If this is not the case (i.e., 
essentially the protostellar luminosity is too weak to modify the 
Jeans length and mass) the relevant length scales and masses for 
collapse revert to the usual equations [Tj and [2] This may occur for 
very low densities and/or high temperature clouds. We have also 
assumed that the mean density of the cloud is low enough that the 
cloud is optically thin on length scales ^> A e a to the infrared ra- 
diation from the protostar. If the density is high enough that the 
radiative feedback is diffusive then using equations [4] and [5] but 
replacing equation [3] with the diffusion approximation 
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the effective Jeans mass can be shown to scale as 
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The interesting thing about this equation for the effective Jeans 



Thus, its dependence on density is somewhat stronger (though still 
weaker than the standard Jeans mass) but the luminosity depen- 
dence is similar. We do not discuss this case further here because 
it applies to densities and/or temperatures that are very different to 
local star-forming regions (e.g. mean densities ^ 10 8 cm -3 and 
temperatures 30 K). 

In summary, it appears from both the numerical simulations 
and the above analytic arguments that radiative feedback may act 
to force a given amount of gas in a dense molecular core to form a 
particular number of stars, with only a very weak dependency of the 
number of stars on the density of the dense core. This process can 
be seen in action in the BBB2003 RT0.5 and BB2005 RT0.5 calcu- 
lations where the main dense cores with masses « 10 Mq in both 
simulations each produce w 10 objects. The characteristic length 
scale between protostars, A c ff , decreases with increasing density in 
such a way as to leave the characteristic mass only weakly depen- 
dent on the cloud's initial density. We propose that this effective 
Jeans mass is responsible for setting the characteristic mass of the 
IMF, and that this may explain the apparent universality of the IMF 
observed in the local Universe. 

Needless to say, the above analytic explanation of how radia- 
tive feedback diminishes the dependence of the characteristic mass 
of the IMF on the initial Jeans mass in the molecular cloud is a 
gross oversimplification. In the calculations, and in reality, molec- 
ular clouds have highly inhomogenous densities, the protostellar 
accretion rates are variable both spatially and with time, and a pro- 
tostar's luminosity depends on its accretion rate, mass, radius, and 
on the details of how discs funnel matter onto the protostar. How- 
ever, the above discussion does provide a framework within which 
we can begin to understand the role of radiative feedback in setting 
the characteristic mass. Larger radiation hydrodynamical calcula- 
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tions that produce many more stars and brown dwarfs to provide 
better statistics will be needed to test whether or not the character- 
istic mass of the IMF does display the above weak dependencies on 
the initial conditions. 



of 0.5 and 5.0 AU). Future calculations should attempt to include 
the additional radiative feedback from within the unresolved re- 
gions surrounding the protostars. 



5 CONCLUSIONS 

We have repeated the previous hydrodynamical star cluster forma- 
tion simulations of |Bate et aL]p003| > and |Bate & Bonnell| l [2005| >, 
but using a realistic gas equation of state and radiative transfer 
in the flux-limited diffusion approximation rather than the origi- 
nal barotropic equation of state. We find that radiative feedback, 
even from low-mass protostars, has an enormous impact on the star 
formation process. Our conclusions are as follows. 

(i) Whereas star formation in the barotropic calculations con- 
tinued unabated in each dense core until the simulations were 
stopped, radiative feedback from newly formed protostars strongly 
suppresses the production of new objects in low-mass (w 10 Mq) 
dense molecular cores after roughly one local dynamical time. 

(ii) Radiative feedback inhibits the fragmentation of massive 
circumstellar discs surrounding newly-formed protostars and any 
dense filamentary gas in their vicinity. This effect, together with 
(i), decreases the numbers of objects formed in the calculations by 
a factor of « 4 compared with the barotropic calculations. How- 
ever, it does not stop the formation of binary and multiple systems. 
Even binaries with separations of only a few AU exist at the end 
of the calculations. The components of these systems are typically 
widely separated when they form, but evolve to close systems via a 
combination of dynamical interactions, gas accretion, and interac- 
tions with discs. 

(iii) The decrease in the fragmentation of discs and dense gas 
near existing protostars results in many fewer dynamical ejections 
than in the barotropic calculations. Since dynamical ejections are 
responsible for the formation of brown dwarfs and low-mass stars 
in these calculations, the radiation hydrodynamical simulations 
produce many fewer brown dwarfs than in the barotropic calcu- 
lations. This results in a ratio of stars to brown dwarfs of w 5:1 that 
is in much closer agreement with observations than the barotropic 
simulations which produced more brown dwarfs than stars. 

(iv) Whereas the characteristic stellar mass was found to scale 
linearly with the initial Jeans of the clouds in the barotropic calcu- 
lations, with radiative feedback the two IMFs are indistinguishable. 
We propose that the reason there is little observed variation of the 
IMF in the local Universe is because the star formation process 
self-regulates itself via radiative feedback. Based on the numerical 
results, we present an analytic argument for how a characteristic 
mass based on radiative feedback from low-mass protostars might 
be derived that scales very weakly with the initial conditions in 
molecular clouds. For example, assuming grey radiative transfer, 
we obtain a characteristic mass that scales as M e g oc p -1 ^ 5 L^ w , 
where p is the density of the cloud and L* is the typical protostellar 
luminosity. 

(v) Finally, we note that due to the sink particle approximation 
used in the radiation hydrodynamical calculations presented here, 
the protostellar luminosity is underestimated. The intrinsic proto- 
stellar luminosity and a substantial fraction of the accretion lumi- 
nosity is neglected. Thus, the dramatic effects of radiative transfer 
presented here are actually lower-limits. We investigate the degree 
to which this may affect our results by performing one of the two 
calculations with different sink particle parameters (accretion radii 
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